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We construct quantum field theory in an analogue curved spacetime in Bose-Einstein condensates 
based on tiie Bogoliubov-de Gennes equations, by exactly relating quantum particles in curved 
spacetime with Bogoliubov quasiparticle excitations in Bose-Einstein condensates. Here, we derive 
a simple formula relating the two, which can be used to calculate the particle creation spectrum 
by solving the time-dependent Bogoliubov-de Gennes equations. Using our formulation, we nu- 
merically investigate particle creation in an analogue expanding Universe which can be expressed 
as Bogoliubov quasiparticles in an expanding Bose-Einstein condensate. We obtain its spectrum, 
which follows the thermal Maxwell-Boltzmann distribution, the temperature of which is experimen- 
tally attainable. Our derivation of the analogy is useful for general Bose-Einstein condensates and 
not limited to homogeneous ones, and our simulation is the first example of particle creations by 
solving the Bogoliubov-de Gennes equation in an inhomogeneous condensate. 
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PACS numbers: 03.75.Kk, 03.75.Hh, 05.30.Jp, 04.62. 

I. INTRODUCTION 

Since the establishment of the theory of general relativ- 
ity by A. Einstein, the study of classical gravity has been 
greatly extended. It provides some fascinating research 
subjects, especially cosmology and black holes. These 
two topics are powerful motivations for studying physics 
in curved spacetime. 

The quantum effects of matter fields in curved space- 
time is also important for studying cosmology and black 
hole physics. In standard modern cosmology, our Uni- 
verse began with an accelerated expansion, known as in- 
flation. The quantum effect of matter fields in the in- 
flationary universe is believed to give rise to cosmologi- 
cal structures such as galaxies. However, this is merely 
theoretical scenario lacking any experimental evidence. 
Another example relates to the very mysterious objects 
known as black holes. According to S. W. Hawking's 
prediction black holes are not completely black but 
rather they emit black body radiation through particle 
creation in the curved spacetime. It is very difhcult to 
confirm this surprising phenomenon experimentally due 
to the extreme low intensity of the emitted radiation. 

With such a background, new attempts to discover 
analogous phenomena in fluid systems have appeared. 
In a pioneering study, Unruh found that Hawking radi- 
ation is not restricted to gravity fields but that it can 
also be observed in a perfect fluid with a sonic hori- 
zon The basic concept involves identifying a sound 
wave propagating in a flowing fluid with a scalar field in 
curved spacetime. In this identification, fiuid flow corre- 
sponds to curved spacetime and provides the spacetime 
metric. The sound wave can then be considered as a field 
in curved spacetime. This analogy has been applied not 
only to Hawking radiation but also to other field theo- 
retical phenomena in curved spacetime, including cosmo- 
logical particle creation [1, 0] . 



One of the merits of considering such an analogy lies in 
possibility of observing quantum effects in a laboratory. 
For the purpose of performing a precise investigation 
of quantum effects, we should consider quantum fluids. 
Bose-Einstein condensates (BECs) in trapped cold atoms 
represent one of the best systems in this context d, H, 0] j 
for the following three reasons: (i) a BEG system is quite 
dilute and weakly interacting, so that its microscopic the- 
ory can be easily studied, (ii) many physical parameters 
of BECs are experimentally controllable, and (iii) various 
quantum phenomena of BECs can be directly observed, 
in stark contrast with other quantum fluids, such as su- 
perfluid helium systems. Therefore, BECs are one of the 
most promising quantum fluids for an analogous system. 

In this spacetime analogy, spacetime corresponds to 
the condensate wave function described by the Gross- 
Pitaevskii (CP) equation ^, [§]. Furthermore, quanta 
of the scalar fleld in spacetime correspond to Bogoli- 
ubov quasiparticles, which are low energy excitations on 
condensates and can be described by the Bogoliubov-de 
Gennes (BdG) equations. Bogoliubov quasiparticles not 
only exist in BECs at flnite temperatures but they are 
also spontaneously created through the motion of BECs, 
which is one of the important problems associated with 
defining a BEC and its superfluidity in nonequilibrium 
Bose systems. 

One of the main subjects of the analogy in BEC is 
to verify the Hawking radiation 111 El, El El, 



El, El, 113, El El IM IM 123. Not only the Hawk- 
ing effect, other aspects of black hole physics as quasi- 
normal frequencies or superradiance were also investi- 
gated [H [II m, m, [23, Hi]. Another stream is on 
cosmological models representin g h omogeneous expand- 



ing universe f2§, 's^, '31', '3^.j33l. |34| or inflationary uni- 
verse [34, 35, 36, 37, 38, 3^, iJ^^The analogy in BEC was 
firstly formulated by Garay et al. and explored by 
Barcelo et al. j4Ij|. In their formulation, the analogue ge- 
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ometry is read from a second-order differential equation 
for the perturbed phase of the condensate wave func- 
tion which obeys the GP equation. In order to connect 
quantum physics in curved spacetime to the behavior of 
a reahstic quantum fluid, Leonhardt et al. investi- 
gated the Hawking effect within the Bogohubov theory 
of the elementary excitations in BEG. More detailed cor- 
respondence between the quanta in analogue spacetime 
and the Bogoliubov quasiparticles was discussed by Jain 
et al. [s^ l. giving an analytical expression for particle 
creation spectrum in terms of the Bogoliubov mode func- 
tions in the case of homogeneous BEG. 

In this paper, we formulate the spacetime analogy in 
terms of the BdG equations as an extension of [ij], by 
exactly relating Bogoliubov quasiparticles with quantum 
particles in curved spacetime. It is shown that the num- 
ber operator of quanta in analogue spacetime is different 
from that of Bogoliubov quasiparticles, unless the cor- 
responding field is normalized correctly. With the cor- 
rect normalization, the number of quanta in analogue 
spacetime can be exactly related with that of Bogoli- 
ubov quasiparticles. Furthermore, it is shown that the 
orthonormal relation between the normalized mode func- 
tions exactly corresponds to that of BdG wave functions. 
We also explicitly define Bogoliubov quasiparticle cre- 
ation as a phenomenon caused by a change in the com- 
plete set of BdG wave functions for the motion of a BEG, 
and we define analogue particle creation in terms of Bo- 
goliubov quasiparticle creation. Our formulation based 
on the BdG equations is generally applicable and gives 
a simple formula for particle creation spectrum in terms 
of the Bogoliubov mode functions. In the derivation of 
the formula, we consider the correctly normalized field so 
that the formula is relevant quantitatively. Our deriva- 
tion of the formula is different from that given in [33 |. 
and consequently the formula is generalized to be appli- 
cable in general cases. Furthermore, our formulation is 
independent of whether condensates are homogeneous or 
inhomogeneous, and the obtained formula is valid even 
for inhomogeneous condensates. 

We then perform a numerical simulation for the sim- 
plest quantum effect in analogue curved spacetime: parti- 
cle creation by the analogue expanding Universe as an ex- 
panding BEG. The spectrum obtained is consistent with 
the thermal Maxwell-Boltzmann distribution, the tem- 
perature of which has already been achieved in actual 
experiments. In this simulation, we set the initial state 
for the excitation field to be vacuum just after the con- 
densate starts to expand, and there is no external ma- 
nipulation after the initial time and the particle creation 
occurs spontaneously, which should be contrasted to the 
previous works [1^, H^, [U, IH, H^l- In this sense, the 
system is isolated. The particle creation simulated here 
is not a result of any perturbation but a field theoret- 
ical phenomenon caused by a difference in definition of 
quasiparticles. 

This paper is organized as follows. In section [TTl we 
derive an analogue spacetime metric from the BdG equa- 



tions. In section [TTTl we construct quantum field theory 
on the analogue spacetime based on BdG theory. In sec- 
tion [iVl we define Bogoliubov quasiparticle creation and 
give a formula for analogue particle creation. In section 
rvl we show the results of a numerical simulation giving 
particle creation in an expanding and shrinking BEG in a 
trapping potential. The conclusions are given in section 

m 



II. ANALOGUE SPACETIME IN BEC 

We start with the Hamiltonian describing bosons 
trapped by an external confining potential 14xt (x) given 
by 



iS: + Vext(x)-Ai 

+ ic/oV'^(x)^(x) 



V'(x). (1) 



Here '0^(x) and '!/'(x) are boson field operators in the 
second quantized form, K = —1t?V^ /2m is the kinetic 
energy operator, Uq is the particle interaction strength 
and /i is the chemical potential [i^,!!^ . In the Heisenberg 
picture, the field equation is given by 



ihdtip{x, t) = [K + yoxt(x) - fi 

+ [/oi^t(x,i)^(x,i)]V'(x,t). 



(2) 



The density operator and the current density operator 
are respectively defined by 



n(x, <) = V'' (x, t)'0(x, t). 



j(x,<) 



2mi 



[^t(x,t)VV^(x,i) 



(3) 



(4) 



In a BEG, the field tjj consists of two components, 
-0(x,t) = ^'(x,t) -I- (^(x, where ^'(x,t) = (V'(x,t)) de- 
scribes the condensate [4J| with (O) denoting a grand 
canonical ensemble average of O. The Bogoliubov field 
(/)(x, t) describes fluctuations about the condensate. The 
condensate wave function 'I'(x, t) obeys the following GP 
equation [1, [§], 

ihdt^i^, t) = [K + K,xt(x) - fi] *(x, t) 

+ C/o|*(x,i)|2^'(x,t). (5) 

The condensate wave function ^'(x, t) is rewritten in 
terms of the amplitude, |5'(x, i)| = ■y/no(x, t), and the 
phase, S{x,t), as follows. 



(6) 



Hereafter, the spatial and time dependence of fields, rep- 
resented by (x, t), is implicit for conciseness of notation. 
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The condensate contributions to the density and current 
density are 

J^o^l^P, jo = ?io— V5 = noVo, (7) 
m 

where vq :— h\7S/m is the velocity of the condensate. 
With these quantities, the GP equation can be rewritten 
as two real equations 

hdtS + n^'/^Knl^^ + Vext - M + UoUo + jv^ = 0, (8) 



dtuo + V • (novo) = 0. 



(9) 



The first equation is Bernoulli's equation for the conden- 
sate current and the second equation represents conser- 
vation of the current. 

Density fluctuations and current density fluctuations 
are described by 

p' = ^*<j) + <j)^^ = ^{e-'^<P + e'^(j)^), (10) 

j' = -^[«'*V0 + 0^V* - (V**)</> - (V0t)*l. (11) 
zmi 

Here, we take the linearized approximation, that is, we 
neglect higher order fluctuations, like (/)''(/). The fluctu- 
ation fields (j), <f)^ obey the BdG equations. These fields 
can be U(l) gauge transformed as 



(12) 



The field equations are 

ihdt4> = (i + 2[/ono + hdtS)4> + UonQ4>\ (13) 
-ihdt4i^ = {L* + 2UonQ + hdtS)^'^ + UqUq^, (14) 

where L is defined as 

L:^k + Vext - ^JL, K -.^ e-'^Ke'^. (15) 
Performing a short calculation, one obtains 

k = K - ihvo - ^ih{W ■ Vq) + ^mvl. (16) 

With these fields, p' and j' take the simple forms 

p' = V^(0 + 4>^), j' = p'vo + nov', (17) 
where we have introduced v' — V<i>' with 

h 



(0-(/.t). 



(18) 



Note that the field $' is the velocity potential in the first 
order approximation with respect to fluctuation fields. 

Taking the difference of eqs. and (fTi|) . we see that 
the operators p' and j' obey the conservation law 



Therefore, the continuity equation is satisfied within first 
order terms for fluctuation fields. Under the hydrody- 
namic approximation, in which the density gradient is 
smooth over the local healing length <^ = h/ \/2mUono, 
the sum of eqs. ((T3)) and fM]) yields 

— (at+vo-V)$'. (20) 
Then, eq. can be rewritten as 

[dt + V- vo) £(9t + vo • V)$' - V • (£c2vl>'^ ~ 0, 

(21) 



where Cs := yX^ono/m. Now, we introduce a symmetric 
tensor 



Un 



-vob 



(22) 



where A is a constant and a,b = 1,2,3. We can then 
write eq. (|2T|) as 



0. 



(23) 



If we identify g^^, with a spacetime metric through ds'^ = 
gf^i^dx'^dx'^ , the field $' obeys the equation of motion for 
a minimally coupled scalar field in curved spacetime with 
the metric (22), and can be seen to be a field in curved 
spacetime. 

The equation for the field $' given in ((2T|) was already 
derived as a perturbation equation of the GP equation, 
see for example |32|, 1451]. In this section, we have red- 
erived it directly from time dependent BdG equations as 
an extension of analysis given in Ts*]. This derivation 
is physically equivalent to the original derivation of ana- 
logue metric given in [lo| . 



III. CONSTRUCTION OF QUANTUM FIELD 
THEORY IN CURVED SPACETIME 

Now, the field $' is found to be a corresponding field 
in analogue spacetime. In this section, however, we show 
that the number of quanta described by $' in the ana- 
logue spacetime does not quantitatively correspond to 
that of quasiparticles in BEG, and there is an exact ex- 
pression for a quantum field in analogue spacetime which 
is different from To address this point, we construct 
quantum field theory on the analogue spacetime by using 
the BdG theory. 

For that purpose, we introduce linear transformation 
and separate operators and wave functions so that the 
operators aj and respectively become annihilation and 
creation operators, as 



V + V-j' = 0. 



(19) 



(24) 
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The operators aj and Oj are set to satisfy the boson 

commutation relation [aj, a^] = Sjk- The wave functions 
Uj and Vj obey the foUowing matrix equations: 



W -Uono 
Uon„ -W* 



(25) 



where W := L + 2C/on-o + tidtS. They then satisfy the 
orthonormal conditions 



d^x {u*Uk - v*Vk) = Sjk, 

(fix {UjUk — VjVk) = 0. 



(26) 



The operators aj and describe the annihilation and 
creation of Bogoliubov quasiparticles, respectively. Here, 
we identify these with annihilation and creation opera- 
tors of quanta in the analogue spacetime. 

The field $' is expanded in terms of operators by use 
of the expansion ([24]) . 



where fj is defined as 



(27) 



(28) 



Firstly, we check that the coefficient functions fj satisfies 
the relativistic wave equation with the metric (|22p . Eq. 
dS]) leads to 



W = K' + UqUq - ihDy, 



with these definitions K' := K ~ tlq [Kitlq^) and 

1 



(29) 



D„ :=vo- V + -(V-vo) 



(30) 



The matrix equation (|25p then gives equations for the 



fimctions u. 



Vj and Uj 



v-i as 



ih [dt + Dy) {uj + Vj) w 2UonQ{uj — Vj), 
ih {dt + Dy) [uj ~ Vj) — K'{uj + Vj), 

where the hydrodynamic approximation, 

K' + 2C/o?io ~ 2Uono, 

has been used. From eq. ([?T|) . 

ih 



(dt + Dy) (uj +Vj). 



' ' 2Uono 
Substituting this into eg. ([5^ . we obtain 



(31) 
(32) 



(33) 



(34) 



D„)^J—{dt + Dy){uj + Vj) + ^K'{uj + Vj) = 0. 



(35) 



The conservation of the condensate current ([9]) leads to 
X 1 



{dt + vo • V) 



{dt + Dy) 



1 



X = 



1 



{dt+Dy)X, 



(at + v-vo)x. 



/no yfno 
for any function X . Therefore, we obtain 

(5t + V-vo)-^(at+vo- V)/,- 

+ ^^K'^}, = 0. 
Furthermore, by noting that 

1 



(36) 
(37) 



V • (noVX), 



(38) 



(39) 



(40) 



with metric (|22|). as expected. 

Now, we calculate the Klein-Gordon (KG) inner prod- 
uct between functions {fj} in the analogue spacetime, 
which is a conserved product for solutions to the wave 
equation in curved spacetime. The KG product for func- 
tions F and G is defined as 



2m y^Jio 

one finds that the functions {/j} satisfy 
^^dcigg'^-d^f,] - 0, 



{F,G) 



KG 



\ [ {Fn>'d^,G* ~ G*n^'df,F)Vd^d^x, 

(41) 



where n — n^d^ is a unit normal vector orthogonal to the 
spacelike surface E, that is, rt is a future-directed timelike 
vector, y/deth d^x is an invariant volume element of the 
surface E. hab is the induced metric on E. If F and G 
are solutions of the Klein-Gordon equation (|40| . the value 
{F, G)kg does not depend on the choice of E as long as S 
is a Cauchy surface 46] . If we choose a i-constant surface 
as E, the induced metric is 



dsl = habdx°'dx^ — CsSabdx"" dx^ 



mA 
'U'o 



(42) 



Then, we obtain Vdeth = {mAcs)^^^!/^ The normal 
vector to the surface is 



^^9, = c;3/^/i^(a, + vo.V). 

" mA 



(43) 



By using eqs. ((36)) and ((3T|) . the KG product between 
the coefficient functions can be calculated as 



ifjJ~> 



Ah 



j,Jk)KG 



d^x [u,ul - v,vl] . (44) 



Using the orthonormal relation (|26p . we obtain 

{fj,fk)KG = Sjk- 

m 



(45) 
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Ah 



Sih and 



Similar calculations show that {f* , f^)KG = —^Ojk 
{fj, fl)KG = 0. Therefore, it is seen that the coefficient 
functions are normal each other. At the same time, it is 
also seen that they are not orthonormal. We find that a 
set of functions 



(46) 



is a set of orthonormal mode functions for the field 
Then the number operator of quanta described by <t' is 



m ~ 



Ah 



Nj — — ala^. 
m •> 



(47) 



The number of quanta in analogue spacetime does not 
correspond to that of Bogoliubov quasiparticles. In order 
to identify these numbers, the field has to be normalized 
so as to make the coefficient functions orthonormal. It is 
easily seen that a set of functions 



(48) 



is a complete set in quantum field theory in analogue 
curved spacetime, that is, these functions satisfy 



{f],fk)KG = Sjk- 

Then the normalized field is 



J2 



Ah 



Then, the number operator of the field ip is 



(49) 



(50) 



(51) 



We have obtained orthonormal mode functions (|48p and 
normalized scalar field (j50p in the analogue spacetime 
from BdG wave functions. Equivalently, we have exactly 
related one Bogoliubov quasiparticle with one quantum 
in the analogue spacetime. Now, with this normalized 
field quantum field theory in curved spacetime can be 
applicable to Bogoliubov quasiparticles quantitatively. 



IV. PARTICLE CREATION 

In this section, we discuss a simple example of particle 
creation and derive a formula connecting particle creation 
spectrum and BdG wave functions. 

Let us consider a condensate which is static and hav- 
ing no excitation at initial time t — t\. For t\ < t < t2, 
the condensate evolves dynamically and, at i = t2, it 
becomes static again. Then, the corresponding effective 
spacetime evolves dynamically and, in general, quantum 
field theory in curved spacetime will predict particle cre- 
ation. Therefore, the analogy implies that, Bogoliubov 
quasiparticles will be created at t = t2. 



We choose a complete set at time t = t\ (A = 1,2), 



{fj^^}, as each /j 



(A) 



the BdG equations uj^'' and Vj^^ at each time t = t\ 



consists of energy eigenfunctions of 

In 

this way, we define one particle state of the quantum by 
using annihilation and creation operators aj^^ and ctj'^^^ 
which diagonalize the BdG Hamiltonian. 

At the initial and final times t = t\ {\ = 1,2), we 
expand the field using a complete set at each time as 



JA)t „(A)* 



(52) 



These two sets are related to each other by linear trans- 
formations 



(53) 



Therefore, the creation and annihilation operators are 
related as 



„(2) _ 



(54) 



Since the initial condensate has no excitations, the state 
for Lp is in the initial vacuum denoted by |0)(i), i.e., 

Q!^^^|0)(i) = 0. We consider the Heisenberg picture and 
the state for ip is not time-dependent. Then, the expecta- 
tion value of the final number operator in the j-th mode 



N. 



(2) 



aj^^^ckj-^'' is calculated as 



(1) 



l^f |0)(i) 



k 



\B 



•]k\ 



(55) 



If Bi^j is not zero for any fc, then particles will be created. 
The expected total number of particles is thus ti{B*B). 

We should note that solutions to the BdG equations uj , 
Vj solve relativistic wave equations (1^0]) through the re- 
lation Therefore, if the time-dependent BdG equa- 
tions are solved, then the functions {/j^''} satisfy the 
equation of motion in time-dependent effective space- 
time. Therefore, time evolution of the initial complete 
set in analogue spacetime is solved automatically. Then, 
the coefficient Bjk is calculated as the KG product at the 
final time t = t2, 



Bjk - - 



d^x{u 



,(2)„(1) 



(56) 



In the derivation of this formula, we have used ((36|) and 
(I5T|) as in the calculation of (^5]) . Here, it should be noted 
that, in the calculation of Bjk, not the set of functions 
{fj} but {fj} has to be considered so that the Bogoliubov 
coefficient Bjk is relevant as a particle creation spectrum 
quantitatively. This expression for the Bogoliubov co- 
efficient (|56p was obtained in [s^ for particle creation 
in homogeneous BEG. Now, the formula has been 
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shown to be vaHd for particle creation in dynamical con- 
densates even when it is inhomogeneous, because, in our 
formulation, any simplification for the Bogoliubov mode 
fimctions is not needed to be assumed. Therefore, our 
formulation is much more powerful. It is powerful enough 
to show that the formula ([56]) is valid for general cases. 
In this way, the spectrum of particle creation (j55p can 
be expressed in terms of the BdG wave functions via the 
formula (|56|) . 

Finally, we comment that the total number of particles 

J d^x (V^^V) = J d^x + {(f,^)), (57) 

is not conserved. dtN 7^ in our model, since the con- 
servation is closed only in the GP equation for the con- 
densate particles; 

dtj d^x\^\^=0. (58) 

To overcome this unphysical result, we have to account 
for the quantum backreaction which is the effect of the 
dynamical development of the Bogoliubov field to the 
BEG. In this paper, we assume that the quantum back- 
reaction can be neglected when the number of the BEG 
particles far exceeds the number of Bogoliubov quasipar- 
ticles. 



V. NUMERICAL SIMULATION 

Now, we show a numerical simulation giving particle 
creation in an expanding and shrinking BEG in the trap- 
ping potential VJ;xt (x) — m{ujlx-^+ujyy'^+ujlz'^)/2, where 
ojx, ujy, and tOz are the trapping frequencies of the poten- 
tial along the x, y, and z-axes respectively. We consider 
the case in which a BEG is tightly trapped only in the 
X-axis; tUx ^ ujy,u!z, and apply the one-dimensional sim- 
ulation by neglecting the spatial derivatives of the y and 
z-axes. We find the stationary solution of the GP equa- 
tion at the trapping frequency Ux = wi, and then change 
the frequency to lUx ~ cuf at t = ti = to make the 
expanding and shrinking BEG. This situation can be in- 
terpreted as a cosmological expansion and contraction, 
because, in the expansion of condensate for example, the 
sound velocity becomes small and the time interval At 
during which a particle on the condensate can travel over 
it. 

At -.^2 / c^^dx, (59) 
Jo 

becomes long, where Rtf is the Thomas- Fermi edge of 
the condensate. This implies that the size of the ana- 
logue spacetime measured by the sound velocity also ex- 
pands. Furthermore, the sound velocity is not constant 
with respect to the spatial coordinate, which implies that 
the analogue spacetime is inhomogeneous. Therefore, 
the analogue spacetime corresponds to an inhomogeneous 



expanding universe. Similarly, a shrinking BEG corre- 
sponds to an inhomogeneous contracting universe. 

Although the condensate starts to move due to chang- 
ing the trapping frequency, we are not interested in par- 
ticle creation caused by such an external manipulation 
but in spontaneous one as a quantum field theoretical 
phenomenon. For that purpose, we set vacuum for Bo- 
goliubov quasiparticles at time just after changing the 
frequency, say t = ti +0. This is the initial condition in 
this simulation. After t — ti+0, therefore, the system is 
isolated without external manipulation, and particle cre- 
ation can be regarded as a consequence of time evolving 
analogue spacetime. 

The spectrum of particle creation X]fcl^j*:P '^''^^ 
calculated as follows. Firstly, we calculate the time evo- 
lution of uj^'' and uj^'' with that of the condensate wave 
function by using eq. starting from the quasi- 

steady solution at t = ii + 0: 

Then, dX t — t2, we obtain u'p and v'^^^ after the time 

evolution and and by using eq. (j60p again. Fi- 
nally, we obtain the spectrum of particle creation from 

eq. dsni). 

As the physical parameters, we use m — 1.44 x 10^^^ 
kg, Wi = 150 X 27r Hz, ujf — cjj/v^, and ujy — ujz — — 
50 X 2tt Hz considering a BEG experiment of ^^Rb atoms, 
the total number of which is set to be iV = 2.5 x 10^. As 
the coupling constant of atomic interactions Uq, we use 
s-wave scattering: 

with a scattering length of a = 5.61 nm. As a numerical 
method, we use the fully dealiased Ghebyshev-Galerkin 
method in space ^Idux < x < 16ax with the Dirichlct 
boundary condition and the Runge-Kutta-Gill method in 
time with At = 1 x lO^^/cji, where Qx = yhlmw-^ is the 
characteristic length defined by the trapping potential. 

Figure [Ha) shows the time development of the 
Thomas-Fermi edge _Rtf of the condensate fitted by a 
harmonic function. Dynamics of the condensate is almost 
periodic, and its period is T ~ 5.10/a;i which is consis- 
tent with that of the breathing mode ^ 2-k j {^^f^j2L0\) . 
The sound velocity Cs decreases with an increase in i?TF 
as shown in Fig. [ijb), also becoming periodic in time. 

The decrease in the sound velocity shown in Fig. [Hb) 
corresponds to expansion of the analogue spacetime. We 
select the time at which the size of condensate becomes 
maximum as the final time <2j i-e., t2 '■— T/2 ~ 2.55/^1, 
and calculate u'^'^ and v^^^ (A = 1,2). Before calcu- 
lating the spectrum of the particle creation, we need to 
verify whether the effect of quantum backreaction can 
be neglected, as discussed above. Figures Ufa) and (b) 
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FIG. 1: Dependence of the Thomas-Fermi edge _Rtf (a) and 
the sound velocity Cs ~ ^/Uono{x — 0)/m (b) at the center 
on time. Here, = ^/W/rnuTi, — ^huji/m are the charac- 
teristic length and velocity defined by the trapping potential. 



show the time dependence of the total number of non- 
condensate particles; 



,(1)|2 



(62) 



and its spectrum at i = t2. Both figures reveal that the 
number of Bogoliubov quasiparticles is small enough rel- 
ative the number of BEC A^o '■— J dx\'^\'^ that the quan- 
tum backreaction will not alter the result qualitatively. 



0.25 r 

0.2 
^=0.15 
0.1 
0.05 

■ 



(a) 



4r 
Cif 2 

S 1 



(b) 







4 6 
t u i 



10 



-20 20 40 60 80 100 



FIG. 2: Dependence of the number of the Bogoliubov quasi- 
particles on time (a) and its spectrum (b). 

Next, we investigate its particle creation. Figure [3] 
shows the spectrum of the particle creation J2k\^jk\'^- 
In the low energy region Ej < /j, corresponding to Fig. 
[Sja), the spectrum has the thermal Maxwell-Boltzmann 

distribution: J2k exp(— i^j^'/ZcBTc) with the 

temperature fc^Te ~ O.lllhiUi. This result means that 
the spectrum of the particle creation in the analogue 
expanding universe obeys the thermal distribution with 
some temperature Te. It is widely believed that the ther- 
mal spectrum of particle creation is caused by formation 
of a static horizon. However, at the time t2, there is no 
sonic horizon. This particle creation is quite nontrivial. 



To confirm this experimentally, a BEC is needed at 
a temperature lower than Tg , to ensure that it is not 
be submerged by thermal noise. Using realistic param- 
eters. To is estimated to be 5.60 nK which has been al- 
ready achieved in many BEC experiments. On the other 
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FIG. 3: The spectrum of particle creation at t = 2.55/a;i in 
the low energy region 5hiOi < Sj'^' < 2QhuJi (a) and lOhuji < 
E'-p < 200huj. (b). 



(2) 

hand, in the high energy region E'j ^ /i shown in 
Fig. [3] (b), the spectrum exhibits power-law behavior 



(2)-4 



This power-law behavior, however. 



seems to be difficult to detect because particle creation 
in this region is very weak. In the high energy region, 
the modes are not hydrodynamic and the analogy breaks 
down. The power-law behavior might be explained by 
trans-Planckian physics in analogv[45i l47j. 

Note that the particle creation shown in this simulation 
occurs spontaneously in the sense that the initial state 
for the excitation field has no excitation at the time just 
after the changing of the trapping frequency, and there 
is no external manipulation after the condensate starts 
to move. This is contrasted to the previous simulation 
giving homogeneous cosmological analogue [3^ , in which 
the external manipulation like changing strength of the 
atomic interaction is required. 



VI. SUMMARY 

In conclusion, we have formulated quantum field the- 
ory in analogue spacetime based on the BdG equations. 
In this analogy, quanta in curved spacetimes are explic- 
itly exactly related with Bogoliubov quasiparticles on 
condensates. It has been demonstrated that, the or- 
thonormal relations for Bogoliubov modes correspond to 
that for quanta in effective spacetime with respect to 
the KG product. We derive a simple formula for the 
particle creation spectrum in terms of BdG wave func- 
tions, which can be applied to simple dynamical evo- 
lution whose initial and final condensations are quasi- 
static. Furthermore, we calculate the particle creation 
in the analogue expanding universe by numerically solv- 
ing time-dependent BdG equations for an expanding 
BEC. The spectrum obtained is consistent with the ther- 
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mal Maxwell-Boltzmann distribution for the temperature 
Tg ~ 5.60 nK, which is experimentally accessible. This 
supports the experimental testing of particle creation in 
an expanding universe by using a BEC. 

Because we neglected the effect of the quantum back- 
reaction from the Bogoliubov field to the BEC, we are 
unable to calculate the precise amount of particle cre- 
ation. To overcome this deficiency, we are extending the 
model to incorporate the quantum backreaction and will 
report on this in near future. Furthermore, numerical 
simulations of particle creations in dynamical BEC corre- 



sponding to Hawking radiation will be reported in future 
studies. 
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